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Abstract 

We present a set of general results on structural features of the g-state Potts model 
partition function Z{G, q, v) for arbitrary q and temperature Boltzmann variable v for various 
lattice strips of arbitrarily great width Ly vertices and length vertices, including (i) cyclic 
and Mobius strips of the square and triangular lattice, and (ii) self-dual cyclic strips of the 
square lattice. We also present an exact solution for the chromatic polynomial for the cyclic 
and Mobius strips of the square lattice with width Ly = 5 (the greatest width for which 
an exact solution has been obtained so far for these families). In the oo limit, we 

calculate the ground-state degeneracy per site, W{q) and determine the boundary B across 
which W{q) is singular in the complex q plane. 
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1 Introduction 



The g-state Potts model has served as a valuable model for the study of phase transitions 
and critical phenomena |jl|, ^. On a lattice, or, more generally, on a (connected) graph G, 
at temperature T, this model is defined by the partition function 

Z(G,g,t;) = J]e-^^ (1.1) 

{cr„} 

with the (zero-field) Hamiltonian 

where ai = 1, ...,q are the spin variables on each vertex (site) i E G; f3 = [ksT)'^] and {ij) 
denotes pairs of adjacent vertices. The graph G = G{V, E) is defined by its vertex set V and 
its edge set E\ we denote the number of vertices of G as = n{G) = \V\ and the number of 
edges of G as e{G) = \E\. We use the notation 

K = pj , a = e^ , v = a-l (1.3) 

so that the physical ranges are (i) a > 1, i.e., f > corresponding to oo > T > for the 
Potts ferromagnet, and (ii) < a < 1, i.e., —1 < f < 0, corresponding to < T < oo for 
the Potts antiferromagnet. One defines the (reduced) free energy per site / = —f3F, where 
F is the actual free energy, via 

f{{G},q,v)= limln[Z(G,g,t;)i/T (1.4) 

where we use the symbol {G} to denote lim„_>oo G for a given family of graphs. In the 
present context, this n — oo limit corresponds to the limit of infinite length for a strip graph 
of the square lattice of fixed width and some prescribed boundary conditions. 

Let G' = {V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex 
set V and an edge set E' C E. Then Z{G, q, v) can be written as the sum 0] 

Z{G,q,v)= J2 q'^'^'^v^'"'^ (1.5) 

G'CG 

where k{G') denotes the number of connected components of G' . Since we only consider 
connected graphs G, we have k{G) = 1. The formula (|1.5| ) enables one to generalize q from 
Z+ to M+ (keeping v in its physical range). The formula ( |1.5|) shows that Z{G,q,v) is a 
polynomial in q and v. 
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The Potts model partition function is, up to a prefactor, equal to a quantity of major 
current mathematical interest, the Tutte (also called Tutte/ Whitney) polynomial. The Tutte 
polynomial of a graph G, T{G,x,y), is given by 

T(G,x,y)= 5^(x-l)'=(^')-^(^)(y-ir(«') (1.6) 

G'CG 

where k{G'), e{G'), and n{G') = n{G) denote the number of components, edges, and vertices 
of G', and c{G') = e{G') + k{G') — n{G') is the number of independent circuits in G'. As 
stated in the text, k{G) = 1 for the graphs of interest here. Now let 

x=l + ^ (1.7) 

V 

and 

y = a = v + l (1.8) 

so that 

g = (x-l)(y-l) . (1.9) 

Then 

Z{G, g, v) = {x- l)'=(^)(y - 1)"(^)T(G, x, y) . (1.10) 
For a planar graph G the Tutte polynomial satisfies the duality relation 

T(G,x,y) = T(G*,y,x) (1.11) 

where G* is the (planar) dual to G. Some reviews of Tutte polynomials include 0-0. 

One interesting special case of the Potts model is the zero-temperature limit (w = — 1) 
of the Potts antiferromagnet, where the partition function is identical to the chromatic 
polynomial P(G, q) counting the number of ways of coloring the vertices of a graph with q 
colors subject to the condition that no adjacent vertices have the same color |1T0| , ^ [11], 12 



Z(G,g,-l) = P(G,g) . (1.12) 

The minimum number of colors necessary for this coloring is the chromatic number of G, de- 
noted x{G)- For sufficiently large g, on a given lattice or graph G, the Potts antiferromagnet 
exhibits nonzero ground state entropy (without frustration). This is of interest as an excep- 
tion to the third law of thermodynamics ||TB| , |T1|. A physical example of residual entropy at 
low temperatures is provided by ice |T^. Nonzero ground state entropy is equivalent to a 
ground state degeneracy per site (vertex), W > 1, since 5*0 = /ebIuVT. We have 

W{{G},q) = lim P(G,g)^/" . (1.13) 
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A subtlety in the definition ( |1.13|) due to the noncommutativity 



hm hm P(G,g)i/" ^ hm hm P(G,g)i/" (1.14) 

g— >(7s n— >oo n— >oo <?— >gs 



at certain points (typically, g = 0, 1, .., x{G)) was discussed in |T6 



Using the formula ( |1.5| ) for Z{G, q, v), one can generalize q from Z_(_ not just to but 
to C and a from its physical ferromagnetic and antiferromagnetic ranges 1 < a < oo and 
< a < 1 to a G C. In particular, we shall be interested here in the case of the T = 
antiferromagnet and the corresponding zeros of Z{G, q, —1) = P{G, q) (called chromatic 
zeros). A subset of these zeros can form a continuous accumulation set in the n — oo limit, 
denoted B. This locus occurs where there is a non-analytic switching between different Xz,G,j 
of maximal magnitude and is thus determined as the solution to the equation of degeneracy 



in magnitude of these maximal or dominant Xz,g/^- As discussed earlier |T7| , pTq , |T9|, ^ 



in the infinite-length limit where the locus B is defined, for a given width and transverse 
boundary condition (free or periodic) B depends on the longitudinal boundary conditions 
but is independent of whether they involve orientation reversal or not. Thus, in the present 
context, for a given Ly, B is the same for the cyclic and Mobius strips. Following the notation 
in |]16i, we denote the maximal region in the complex q plane to which one can analytically 



continue the function W{{G},q) from physical values where there is nonzero ground state 
entropy as Ri . The maximal value of q where Bq intersects the (positive) real axis is denoted 
1c{{G}). Thus, region Ri includes the positive real axis for q > qc{{G}). 

In the present work we shall give a number of general results on structural features of the 
g-state Potts model partition function Z{G, q, v) for arbitrary q and temperature Boltzmann 
variable v for various lattice strips of arbitrarily great width Ly vertices and length 
vertices. We shall also discuss an exact solution for the chromatic polynomial for the cyclic 
and Mobius strips of the square lattice with width Ly = 5 (the greatest width for which 
an exact solution has been obtained so far for these families). In the L^ oo limit, we 
calculate the ground-state degeneracy per site, W{q) and determine the boundary B across 
which W{q) is singular in the complex q plane. 

There are several motivations for this work. We have mentioned the basic importance of 
nonzero ground state entropy in statistical mechanics. From the point of view of rigorous 
statistical mechanics, exact solutions are always valuable for the insight that they give into 
the behavior of the given system under study. The results derived here give further insight 
into the general structure of the Potts model partition function for lattice strips. The study of 
VrdG}, q) with q generalized to complex values enables one to gain a deeper understanding 
of the behavior of this function for q G Z+, in much the same way as the study of functions 
of a complex variable imparts a deeper understanding of functions of a real variable in 
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mathematics. In particular, one sees how the value qc on the real axis corresponds to an 
intersection of the locus B with this axis. In addition to the papers [0, 0, |l^, |l^, 
some relevant studies of Potts model partition functions and/or chromatic polynomials for 



strip graphs of regular lattices include |2T|-[^. Related mathematical papers, in addition 



to those already mentioned, include [^-|58[, and Ref. contains further references. 



2 Structural Results 



2.1 General 

In this section we briefly review some general structural results which will be used for our 
new results below. A general form for the Potts model partition function for the strip graphs 
Gm considered here, or more generally, for recursively defined families of graphs comprised 
of m repeated subunits (e.g. the columns of squares of height Ly vertices that are repeated 
Lx = m times to form an x Ly strip of a regular lattice with some specified boundary 



conditions), is |]T9| 



Z{Gm,q,v)= Cz,G,j{Xz,G,j{q,v)y 



(2.1.1) 



where Nz,g,x, ^z,G,j, and Xz,G,j depend on the type of recursive family G (lattice structure 
and boundary conditions) but not on its length m. For strips with periodic longitudinal 
boundary conditions (e.g., cyclic and torus strips) eq. ( p.l.lD follows from the property that 
Z{Gm, q, v) can be expressed as the trace of a certain transfer matrix. 



Z{Gm, q, v) = Tr(T^ 



(2.1.2) 



Let 



Z,G,X 

Cz,G = ^ Cz,G,j ■ 



(2.1.3) 



For the strips with periodic longitudinal boundary conditions, Cz.g is thus the dimension of 
the matrix Tz- In particular, for the cyclic strip graphs of the square and triangular lattices 
considered here, the repeated subgraph is the path graph consisting of Ly vertices, and so, 
denoting Gz,g = Cz,Ly as before, one has |^ 



Cz,L, = q'^ 



(2.1.4) 



For strips with reversed-orientation periodic longitudinal boundary conditions (Mobius and 
Klein bottle boundary conditions), (|2.1.2| ) is modified by the insertion of the requisite 
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orientation-reversal operator in the trace, with the result that |^| 



C 



Z,Mh,Ly — Q 



(2.1.5) 



where [u] denotes the integral part of u. 

Equivalently, for the Tutte polynomial of a recursive graph 



T{Gm,x,y) 



X 



where 



(2.1.6) 



(2.1.7) 



For the special case of the T = antiferromagnet, the partition function, or equivalently, 
the chromatic polynomial P{Gm,, q) has the corresponding form ||23[ 



Let 



Np.G.X 

P{G^,q)= Cp^G,A^p,GM)r ■ 

Np.G.X 

Cp^G = ^ Cp^Gj • 



(2.1.8) 



(2.1.9) 



Here, the strip graph G is indexed by its width, Ly. The analogue of ( p.l.4|) for the chromatic 
polynomial for cyclic strip graphs of the square and triangular lattices is 

Cp,L, = P{Tl,, q) = q{q - if^'' (2.1.10) 
where T„ denotes the tree graph with n indices, and the analogue of (|2.1.5| ) for the chromatic 



polynomial of the Mobius strip of the square lattice is 

for even L,, 



P,Mb,L 



y 



P{T^ Ly+\ y q) for odd Lj, 



(2.1.11) 



2.2 Cyclic and Mobius Strips of the Square and Triangular Lattice 

In it was shown that for cyclic and Mobius strips of the square lattice of fixed width Ly 
and arbitrary length (and also for cyclic strips of the triangular lattice) the coefficients 
Cj{q) in the Potts model partition function (|2.1.1|) , and hence, a fortiori, in the chromatic 
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polynomial (|2.1.8| ), are polynomials in q with the property that there is a unique polynomial, 
denoted c^*^-*, of degree d in q. Further, this was shown to be a Chebyshev polynomial of the 
second kind: 



d-j 



(2.2.1) 



f/2.(g^/V2) = 5^(-l)M . - )g 

j=0 

where Un{x) is the Chebyshev polynomial of the second kind. A number of properties of 
these coefficients were derived in for our present work, we shall need the following 

special values: 



If q=l then c^'^^ 



If q = 2 then c' 



id) 



for q = 

1 a d = mod 3 

if d = 1 mod 3 

-1 iid = 2 mod 3 

1 ifd = 0,l mod 4 
-1 ifd = 2,3 mod 4 



(2.2.2) 
(2.2.3) 

(2.2.4) 



In contrast, for the Mobius strips of the triangular lattice, the coefficients in the Potts model 
partition function and chromatic polynomial do not, in general, have the form c'-'^-' p6| . 

Thus, the terms Xz,Ly,j (equivalently, XT,Ly,j) and XpXyj that occur in ( p.l.l|) , (|2.1.6| ), 
and ( p.l.8|) can be classified into sets, with the \z,Ly,j{q, v), XT,Ly,j{x, y), and \p,Ly,j{q) in the 
d'th set being defined as those terms with respective coefficients cz,Ly,j = CT,Ly,j and cp^Lyj 
being equal to c^'^\ In Ref. the numbers of such terms, denoted nz{Ly, d) and np{Ly, d) 
respectively, were calculated. Note that 



nz{Ly,d) = nT{Ly,d) . 



(2.2.5) 



Labelling the eigenvalues with coefficient c^'^^ as \z,Ly,d,j with I < j < nz{Ly,d), the Potts 
model partition function can be written in the form 



nz(Ly,d) 



3=1 



Z{G[LyXm,cyc.lq,v) = J2\c^'^ E 
or equivalently, for the Tutte polynomial, 

T{G[Ly X m,cyc.],x,y) 

where 



Z,Ly,d,j) 



(2.2.6) 



X 



d=0 



nT{Ly,a 



{^T,Ly,d,j) 



^Z,LyA,M^^) = V^'>'T,Ly,d,ji^^y) 



(2.2.7) 



(2.2.8) 



with the relations (11.71) and (11.81). 
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With analogous labelling, the chromatic polynomial has the form 



(2.2.9) 



Combining ( p.2.6| ) with ( p.l.4| ) and ( |2.2.9| ) with ( |2.1.10| ), one has the following relations 
for the cyclic strip graphs of the square and triangular lattice of width Ly and arbitrary 
length H 



d=0 



Cp,L, = Y.c^''^np{Ly,d) = q{q-lfy 



(2.2.10) 



(2.2.11) 



d=0 



The total numbers, Nz,Ly,x and Np^^^^x, of different terms Xz,Ly,j and Xp,Ly,j in eqs. ( |2.1.1D 
and ( |2.1.^ ) are given by P5| 

Ly 

Nz,Ly,x = NT,Ly,x = J2 M^y^ d) (2.2.12) 



and 



Np,Ly,x = ^np{Ly,d) . 



(2.2.13) 



d=0 



We shall recall here some results on the numbers np{Ly,d)] these are nonzero for < 
d < Ly and are given by 

np{Ly,Ly) = l (2.2.14) 
np(l,0) = l (2.2.15) 
with all other numbers np{Ly, d) being determined by the two recursion relations 

np{Ly + 1,0) = np{Ly,l) (2.2.16) 

np{Ly + l,d) = np{Ly,d—l)+np{Ly,d)+np{Ly,d+l) for Ly > 1 and 1 < d < Ly + 1 . 

(2.2.17) 

The solution to these relations yielded, in particular, the results 



npiLy,0) = MLy^ 



where M„ is the Motzkin number, 



i=o ^ 



(2.2.18) 



(2.2.19) 
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where 



is the Catalan number. Further, 



and 



n + l\ n 



2n 



np{Ly, 1) = Ml, 



(2.2.20) 
(2.2.21) 



np{Ly,Ly-l) = Ly . (2.2.22) 

The total number of Xz,Ly/s in the Potts model partition function Z([Ly x m, eye], q, v) or 
Z{[Ly X m,Mb],q,v) is ^ 

Nz,L,,x = (^^'] (2.2.23) 



and the total number of Xp^Lyj^^ in the chromatic polynomial P{[Ly x m, eye], q) or P{[Ly 
m,Mh],q) is [||] 



iVp,..,. = 2(L, - 1)! Ey,).,[^i2^.), 



(2.2.24) 



For arbitrary Ly, eq. (|2.2.14|) shows that there is a unique \p,Ly,d corresponding to the 
coefficient e^^^ of highest degree, d = Ly, and this term is [41 



A 



P,Ly,d=Ly,j = l 



(2.2.25) 



As indicated, since this eigenvalue is unique, it is not necessary to append a third index, 
as with the other A's, and we avoid this for simplicity. The result ( p.2.25|) corresponds to 
the analogous property that in the full Potts model partition function or equivalent Tutte 
polynomial, nz{Ly,Ly) = 1 and the unique Xz,Ly,d=Ly with coefficient c'-'^^^ in Z{[Ly x 
L^,eye.],q,v) or Z{[Ly x L^,Mh],q,v) is 



A; 



(2.2.26) 



Equivalently, the unique \T,Ly,d=Ly with reduced coefficient cr,G,i ~ '^''^^^ T{\LyXLx, eye.], x, y) 
or T{\Ly X L^,Mh],x,y) is 

\T,Ly4=Ly = 1 . (2.2.27) 

(Here, one should not confuse the arguments of the Tutte polynomial x, y with the longitu- 
dinal and transverse directions x,y.) 
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For a strip of fixed width Ly and arbitrary length = m with Mobius boundary con- 
ditions, {FBCy,TPBCx), the general form of the Potts model partition function was deter- 
mined in [E6i to be 



ZiG[LyXm,Mb],q,v)= 



nz,Mb{Ly,d,'q) 



d=0 



r)=±l 



"Y i>^z,Ly,d,vjy 



where 



d 



max, Mb 



-Y for even Ly 
for odd Ly 



(2.2.28) 



(2.2.29) 



Correspondingly, 

^max .Alb 

T{G[LyXm,Mb],x,y) = ^ {^^"^ E ^ E (^W,,,,)'"] } (2.2.30) 



nT,Alb{Ly,d,v) 



X 



d=0 ri=±l j=l 

with nT,Mb{Ly, d, rj) = nz,Mb{Ly, d, rj). For the special case of the chromatic polynomial, 

dmax,Mb np,MbiLy,d,'q) 

P{G[LyXm,Mb],q)= E {^^'^ E ^ E (^W,..)"] } • (2-2.31) 



d=0 



ri=±l 



Combining ( [^.2. 28] ) and ( |2.2.31| ) with these formulas, one has the relations for the Mobius 
strip graphs of the square and triangular lattice of width Ly and arbitrary length 

^■max ,AIh 

Cz,Mb,Ly = E ^^'^^ E Vnz,Aib{Ly, d, 7]) = g^^^ 

d=0 r]=±l 



(2.2.32) 



Cp,Mb,Ly = 2^ c^' 2^ r]npMLy,d,r]) = i p{T q) fo^. odd Ly ■ (2-2-33) 

In the nzMb{Ly,d,ri) and np^Mb{Ly,d,ri) were determined and a set of transformations 
were given that relate the Potts model partition functions for the cyclic and Mobius strips. 



2.3 Self-Dual Cyclic Strips of the Square Lattice 

Consider a strip of the square lattice with (i) a fixed transverse width Ly, (ii) an arbitrarily 
great length L^, (iii) periodic longitudinal boundary conditions, such that (iv) each vertex 
on one side of the strip, which we take to be the upper side (with the strip oriented so that 
the longitudinal, x direction is horizontal) are joined by edges to a single external vertex. A 
strip graph of this type will be denoted generically as (where the subscript D refers to the 
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self-duality) and, when its size is indicated, as GoiLy x L^). Although this family of graphs 
differs from the simple recursive cyclic or Mobius lattice strips, owing to the feature that all 
of the vertices on the upper side are connected to a single external vertex, we showed ^0|, ^ 
that the Potts model partition function and its special case, the chromatic polynomial, for 
the GoiLy X Lx) family have a structure analogous to that for the cyclic and Mobius strips 
of the square lattice and cyclic strip of the triangular lattice, in the sense that there is a 
unique coefficient which is a polynomial of degree d in q, denoted k^'^\ There are nz,Go{^y^ ^) 
terms Xz,GD,Ly,d,j ^{GniLy x Lx),q,v) having coefficient The Potts model partition 
function, Tutte polynomial, and chromatic polynomial were shown to have the following 



forms f50, 51 



nz,Gr,iLy,d) 



Z{GD[Ly X L,],q,v) = J2 E (^z^o^M: 



(1=1 



(2.3.1) 



or equivalently, for the Tutte polynomial. 



Ly + l 



nT,Gr,{Ly,d) 



T{GD[Ly Xm],X,y)=J2 ['^^'^ J2 (^T,Gn,Ly,cl. 

d=l j=l 



(2.3.2) 



where 



and 



nT,GDiLy,d) = nz,GD{Ly,d) 



(2.3.3) 



K 



id) 



\fq U2d- 

d-l 
j=0 



2d- 



1-J 



3 



-,d-j 



(2.3.4) 



where Un{z) is the Chebyshev polynomial of the second kind. Thus, = q, 



K 



(2) 



g(g-2). 



^^(3) = q(^q _ x)(^g _ and so forth for higher values of d. From ( p..lO|) it follows that 



Xz,GD,Ly,dM^^) = V^'^T,GD,LyAji^^y) 



l<d< Ly + l 



(2.3.5) 



with the relations (|1.7|) and ( |1.8|) . Clearly, for the total numbers, Nz,GD,Ly, 
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For the chromatic polynomial, 



Ly + l 



P{Gn[Ly X LJ, q)=Y, E ^^^p^g^M^Y 



d=l 



(2.3.6) 



The analogues to ( p.2.10|) and ( p.2.11|) for this self-dual strip of the square lattice are 

(2.3.7) 



Ly + l 

E K,'-'^^nz,GniLy,d) = q 

d=l 



Ly + l 



Ly + l 



E K^^^np,G,{Ly,d) = P{TL^+i,q) = g(g - 1)^^ • 



(2.3.8) 



d=l 



A number of properties of the k*^*^^ coefficients were derived in |^|; for our present work, 
we use the following special values 



^'^^ =0 for q = 



iid = mod 3 
If q = l then k^'^^ = { 1 ii d = 1 mod 3 



If q = 2 then k^'^'> = | 



-1 ifrf = 2 mod 3 

if d is even 
l'' if d is odd and d = 2k + 1 



(2.3.9) 
(2.3.10) 

(2.3.11) 



The numbers nz^Goi^y^ d) and np^Goi^y^ d) were determined in [5^] and related to nz{Ly, d) 
and np{Ly, d) for cyclic strips of the square and triangular lattices. For the total number of 
terms in the Potts model partition function. 



^2Ly + r 

The analogous total number of terms in the chromatic polynomial is given by 

1 + X\V2 



(2.3.12) 



l-3a; 



1 



X 



E Np^Go^LyAx"^'^' 



(2.3.13) 



and satisfies 



Np,Gn,Ly,X — -Np^Ly+l,X 



(2.3.14) 



where the latter refers to cyclic or Mobius strips of the square or triangular lattice. 
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3 Structural Sum Rules 



3.1 Relations Between Sums of nz{Ly,d) 

In this section we derive new structural relations for the Potts model partition function 
and chromatic polynomial. These involve sums of terms and hence, following a traditional 
nomenclature in physics, we shall refer to them as sum rules. 

Proposition 1 For the strip of the square or triangular lattice of width Ly and arbi- 
trary length with {FBCy, PBC^) (i.e., cyclic) boundary conditions, 

J2 nz{Ly,d) = J2 MLy,d) = U^^') . (3.1.1) 

0<d<Ly 0<d<Ly ^ y ^ 

d even d odd 



Proof Evaluating ( 2.2.101 ) for g = and using ( p. 2. 2 ) yields the first equality. From 



this equality and the fact that the sum of the number of eigenvalues Xz,Ly,d,j with coefficients 
of even and odd d is equal to (|2.2.23| ), the last equality in (|3.1.1|) follows. □ 

Clearly, eq. ( |3.1.1|) also applies with nz{Ly, d) replaced by riT^Ly, d). Note that the right- 
hand side of the last equality in (|3.1.1D is equal to the number of directed animals of length 



n = Ly on the triangular lattice The first few values of this quantity for 1 < Ly < 5 
are 1, 3, 10, 35, and 126, respectively. 

It is of interest to observe a related structural feature for cyclic strips of the square or 
triangular lattice. We start by noting that for any graph G, Z{G,q,v) = if g = 0, as 
is clear from ( |1.5| ). We apply this to the case where G is the cyclic strip of the square or 
triangular lattice, denoted for short simply as [Ly x Lx,cyc.] with L^ = m, and use ( |2.2.2D 
to obtain 

Z{[Ly X m, eye], q = 0,v) = 

nz{Ly,d} nz{Ly,d) 

= E E (Aw.(o,^)r- E E {^z,Ly,d,jMr . (3.1.2) 

0<d<Ly j=l 0<d<Ly j=l 

d even d odd 

We observe from our explicit calculations that the terms Xz,Ly,d,j{q,v) do not, in general, 
vanish at g = 0. Since the cancellation between the various m'th powers of the eigenvalues 
must occur for arbitrary m, they must, in fact, cancel in pairs. Note that this provides a 
second proof of the result ( p.l.lj) , since the cancellation requires that the number of eigen- 
values Xz,Ly,d,j with even d must be equal to the number of eigenvalues Xz,Ly,d,j with odd 
d. 
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Next, we have 

Prop. 2 For the cychc strip of the square or triangular lattice of width Ly and arbitrary 
length Lj: 

nz{Ly,d)- ^ nz{Ly,d) = l. (3.1.3) 

0<d<Ly 0<d<Ly 

d=0 mod 3 d=2 mod 3 

Proof To prove this, we evaluate (|2.2.1CI| ) for q = I and use (|2.2.3| ). □ 

Prop. 3 For the cyclic strip of the square or triangular lattice of width Ly and arbitrary 
length L3. 

J2 MLy, d)- Yl ^^(^?/' ^) = 2'''' • (3-1-4) 

0<d<Ly 0<d<Ly 

d=0,l mod 4 d=2,3 mod 4 

Proof To prove this, we evaluate (|2. 2.101 ) for q = 2 and use (|2.2.4| ). □ 

Prop. 4 For the cyclic strip of the square or triangular lattice of width Ly and arbitrary 
length L^ 

1 /2L. 



d=0,l mod 4 



0<d<Ly 

d=2,3 mod 4 



(3.1.5) 



Proof This result follows from ( |3.1.4D and the fact that the sum ^ 



jO<d<L, 



nz{Ly,d) 



Nz,Ly,x, for which one has the expression given in eq. ( p.2.23|) . □ 



3.2 Relations Between Sums of np{Ly,d) 

Prop. 5 For the cyclic strip of the square or triangular lattice of width Ly and arbitrary 
length Lj. 

E np{Ly,d) = E npiLy,d) 



0<d<Ly 
d even 



0<d<Ly 

d odd 



{Ly-l)\ E 



[Ly — j] 



^(j!)2(L,-2j)! 



(3.2.1) 
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Proof Evaluating ( 2.2. 11| ) for g = and using ( p. 2. 2 ) yields the first equality. From 
this equality and the fact that the sum of the number of eigenvalues Xp,Ly,d,j with coefficients 
of even and odd d is equal to (|2.2.24| ), the last equality in (|3.2.1|) follows. □ 

Note that the right-hand side of (|3.2.1| ) is equal to the number of directed animals of 
length Ly on the square lattice The ffist few values of this quantity for 1 < Ly < 5 are 
1, 2, 5, 13, and 35, respectively. 

Again, one can remark on a related structural feature for cyclic strips of the square or 
triangular lattice. Clearly, P{G, g) = if g = 0. We apply this to the case where G is the 
cyclic strip of the square or triangular lattice, denoted for short simply as [Ly x Lx,cyc.] 
with Lj. = m, and use ( p.2.2|) to obtain 



P{[Ly X m, eye], q = 0) = 

np{Ly,d) np(Ly,d) 

E E (w.(? = o)r- E E (w.(^ = o)r- (3-2.2) 



0<d<Ly j=l 

d even 



0<d<Ly j=l 
d odd 



We observe from our calculations that the Xp^iy^s are not, in general, zero for q = 0. Since 
the cancellation between the various m'th powers of the eigenvalues must occur for arbitrary 
m, they must, in fact, cancel in pairs. Hence, the number of eigenvalues \p^Ly,d,j with even 
d must be equal to the number of eigenvalues Xp^Ly,d,j with odd d. This yields another way 



of seeing the result (|3.2.1 
Next, 

Prop. 6 For the cyclic strip of the square or triangular lattice of width Ly > 2 and 
arbitrary length L^ 



E np{Ly,d)= E np{Ly,d) 



0<d<Ly 
d=0 mod 3 

Np^LyA — Np^Ly^i^d 



0<d<Ly 

d=2 mod 3 



(3.2.3) 



where the expression for Np^i^^d was given in eq. (|2.2.24|) . 

Proof To prove the first equality, we evaluate ( p.2.11|) for q 



the second equality, we observe that from the results of 

^ np{Ly,d) = Np^Ly-l,\ 



1 and use (^T^. For 
(3.2.4) 



0<d<Ly 

d=l mod 3 
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where the expression for Np^L^^x was given in ( p.2.24|) . Together with the first equahty in 
( p.2.3|) and the formula from |^ for the total sum, (|2.2.24| ), this yields the second equality. 



□ 



For 2 < Ly < Q the values of the right-hand side of the second equality in ( 3.2.3 ) are 1, 
3, 8, 22, 61. The case = 1 is special since for this case Cp^Ly=i = 1 rather than at 
g = 1, so that the analogue to eq. ( |3.2.3|) reads simply rap(l, 0) = 1. 

Prop. 7 For the cyclic strip of the square or triangular lattice of width Ly and arbitrary 
length Lx 

np{Ly,d)- np{Ly,d)=2. (3.2.5) 

0<d<Ly 0<d<Ly 

d=0,l mod 4 d=2,3 mod 4 

Proof To prove this, we evaluate ( 2.2.11 ) for g = 2 and use ( 2.2.4| ). □ 

Prop. 8 For the cyclic strip of the square or triangular lattice of width Ly and arbitrary 
length Lx 



J2 np{Ly,d) = ^Np,Ly,x + l 



0<d<Ly 

d=OA mod 4 



J2 np{Ly,d) = ^Np^L^^x-l. (3.2.6) 



2 

0<d<Lj; 

d=2,3 mod 4 



Proof This result follows from (|3.2.5| ) and the fact that the sum Ylo<d<L ^pi^yi d) 



Np,Ly,x, for which one has the expression given in eq. ( |2.2.24| ). □ 

For an arbitrary connected graph with at least two vertices, P{G,q) = for q = 1. It 
follows, in particular, that 

P{[Ly X m,cyc.],q = 1) = 

np{Ly,d) np{Ly,d) 

= E E i^p,Ly,dA<i = - E E (Ap,L„.„(g = i))"^. (3.2.7) 

0<d<Ly j=l 0<d<Ly j=l 

d=0 mod 3 d=2 mod 3 

However, in contrast to the case for q = 0, the cancellation here is not pairwise because some 
of the A's vanish. Specifically, as will be shown below, the term Xp^Ly,d,j ior d = Ly — 1 and 
j = 1 is g — 1, which vanishes at g = 1. 
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4 Sum Rules Applicable to Mobius Strips of the Square 
Lattice 



4.1 Relations Between Sums of nz,Mb{Ly, d) 

Prop. 9 For the Mobius strip of the square lattice of width Ly and arbitrary length 

[nz,MbiLy, d, +) - nz,Mb{Ly, d, -)] 



d even 



[nz,MbiLy,d,+) - nz,Mb{Ly,d,-)] . (4.1.1) 



0<d<dmax,Mb 

d odd 



Proof To prove this, we evaluate ( 2.2.32| ) for g = and use ( 2.2.2| ). □ 



Applying the general fact that for any graph G, Z{G, q,v) =0 for g = 0, to the Mobius 
strip of the square lattice and using (|2.2.2| ), we obtain 

Z{sq[Mb, LyXm],q = 0,v)=0 

(Az,L,d,,,(g = 0)r- Yl (Az,L,d,,,(g = 0)r 

(4.1.2) 



'JSdSa^ax,Mb "^''■^''■max .Mb 

r]=±l; C=+l r?=±l; C=-l 



where 

C = r]{-iy. (4.1.3) 

Since the cancellation between the various m'th powers of the eigenvalues must occur for 
arbitrary m, they must, in fact, cancel in pairs. This provides another way of understanding 
the result in ( fl.l.lQ . 

Prop. 10 For the Mobius strip of the square lattice of width Ly and arbitrary length 



[nz,MbiLy, d, +) - nz,Mb{Ly, d, -)] 

d=0 mod 3 

+ Y [~^z,Mb{Ly, d, +) + nz,Mb{Ly, d, -)] 



0<<l<<lmax,Mb 

d=2 mod 3 



1 . (4.1.4) 
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Proof To prove this, we evaluate (|2.2.32| ) for q = 1 and use (|2.2.3| ). □ 

Prop. 11 For the Mobius strip of the square lattice of width Ly and arbitrary length 



[nz,Mb{Ly, d, +) - nz,Mb{Ly, d, -)] 

d=0,l mod 4 

+ X] [-nzMbiLy, d, +) + nz,Mb{Ly, d, -)] 

=2,3 rr 

2[^] 



d=2,3 mod 4 



(4.1.5) 



Proof To prove this, we evaluate (|2.2.32| ) for q = 2 and use (|2.2.4| ). □ 



4.2 Relations Between Sums of np^Mb{Ly, d) 

Prop. 12 For the Mobius strip of the square lattice of width Ly and arbitrary length 

[np^Mb{Ly, d, +) - np^Mb{Ly, d, -)] 



d even 



[rap,Mf.(-^^y, c?, +) - np^Mb{Ly, d, -)] . 



(4.2.1) 



d odd 



Proof To prove this, we evaluate (|2.2.33| ) for q = and use ( |2.2.2| ). □ 

Applying the general fact that for any graph G, P{G, g) = for g = 0, to the Mobius 

strip of the square lattice and using ( p.2.2| ), we obtain the following relation, where the A's 

are evaluated at g = 0: 



P{sq[Ly X m,M6],g = 0) = 



E 

0<d<dj-aax,Mb 
rj=±l- C=+l 



(4.2.2) 



0<<l<<lmax,Mb 
ri=±l; C=-l 



where C was defined in ( |4.1.3D . Since the cancellation between the various m'th powers 
of the eigenvalues must occur for arbitrary m, they must, in fact, cancel in pairs. This is 
another way of understanding the result ([4.2. 1|) . 
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Prop. 13 For the Mobius strip of the square lattice of width Ly and arbitrary length 



[np^Mb{Ly, d, +) - np^Mb{Ly, d, -)] 

d=0 mod 3 

['^P.Afbl-^y, d, +) — np^Mb{Ly, d, — )] 



(4.2.3) 



d=2 mod 3 



Proof To prove this, we evaluate ( 2.2.33| ) for q = 1 and use ( 2.2.3 ). □ 

Prop. 14 For the Mobius strip of the square lattice of width Ly and arbitrary length 

[np^Mb{Ly, d, +) - np^MbiLy, d, -)] 

d=0,l mod 4 

+ [-np^MbiLy,d,+) +np^Mb{Ly,d,-)] 



d=2,3 mod 4 



2 . 



(4.2.4) 



Proof To prove this, we evaluate (|2.2.33|) for q = 2 and use (|2.2.4|). □ 



5 Sum Rules for Self-Dual Strips of the Square Lattice 

5.1 Relations Between Sums of nz^Goi^y^^) 

We first note that for the self-dual strips of the square lattice Go^Ly x L^), there is no 
analogue of the sum rules (|3.1.1|) and (|3.2.1|) because all of the coefficients k,^^^ have g as a 
factor and hence vanish at g = 0. We do, however, find several sum rules. 



Prop. 15 For the self-dual strip of the square lattice GniLy x L^] 

nz,GniLy,d)- nz,GniLy,d) = 1 . 



(5.1.1) 



l<d<Ly+l 

d=l mod 3 



l<d<Ly+l 

d=2 mod 3 



Proof To prove this, we evaluate (|2.3.7| ) at g = 1 and use ( p.3.10|) . □ 
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Prop. 16 For the self-dual strip of the square lattice, Go^Ly x L^), 

i-l)'nz,GniLy,d) = 2'^y . (5.1.2) 



l<d<Ly+l 

d=2k+l 



Proof To prove this, we evaluate (|2.3.71 ) at g = 2 and use ( p. 3. 111) . □ 



5.2 Relations Between Sums of np^Goi^yid) 

Prop. 17 For the self-dual strip of the square lattice Go^Ly x L^), 

^ np^GoiLy,d)= np^GniLy,d) . (5.2.1) 



l<d<Ly+l l<d<Ly+l 

d=l mod 3 d=2 mod 3 



Proof To prove this, we evaluate (|2.3.8| ) at g = 1 and use ( p.3.10|) . □ 



Prop. 18 For the self-dual strip of the square lattice Gni^Ly x L^), 

Y i-l)WGniLy^d) = l. (5.2.2) 



l<d<Ly+l 

d=2k+l 



Proof To prove this, we evaluate (|2.3.8| ) at g = 2 and use ( p.3.11|) . □ 



As before for other strip graphs, one may analyze how cancellations occur in the chromatic 
polynomial at values of q where it vanishes. Thus, let us consider the value g = 1 and use 
the fact that for any connected graph G with at least two vertices, P{G, g) = if g = 1. 
Applying this to G = Go^Ly x L^) and using the property ( p.3.10|) , we obtain the following 
relation, where the A's are evaluated at g = 1: 

P{GD{Ly X m,g = 1) = 

np,Gij{Ly,d) np^Gij{Ly,d) 

= X] {^P,GD,LyAj)'^ - X] {^P,GD,Ly,d,j)"^ ■ (5.2.3) 

l<d<Ly+l j=l l<d<Ly+l j=l 

d=l mod 3 d=2 mod 3 

For the cases where we have carried out explicit calculations, we find that the A's do not 
vanish at g = 1 for this family of graphs. Since the cancellation in ( |5.2.3| ) between the 
various m'th powers of the eigenvalues must occur for arbitrary m, they must, in fact, cancel 
in pairs. 
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Since the family GoiLy x L^) contains triangles, P^GoiLy x Lx],q) = for g = 2. One 
may investigate what is implied by this property. Using (|2.3.11| ), one finds 



(5.2.4) 



P{GD[Ly X m],g = 2) = 

np(GD,Ly,d) 

l<d<Ly+l j=l 

d=2k+l 

However, this cancellation does not, in general, occur in a pairwise manner, because some 
of the \p,GD,Ly,2k+i,j vanish. For example, for Ly = 2, np^Gn, 1) = 2 while np{GD,3) = 1, 
but the Xp,Gd,2,i,2, given by eq. (6.3) in [Q, 

1 
2 



g2 - 5g + 7 ± (g^ - 6q^ + ISg^ - 22g + 17) 



1/2 



J = 1,2 



(5.2.5) 



vanishes at g = 2. 



6 Sum Rules for Other Recursive Families of Graphs 

Proceeding in a manner similar to that above, one can obtain analogous sum rules for other 



recursive families of graphs. We comment briefly on one such family. In |5^ we presented 
theorems determining the structure of the Potts model partition function for a cyclic clan 
graph (CG). We recall some relevant definitions. A complete graph Kr is a graph containing 
r vertices with the property that each vertex is connected by edges to every other vertex. 
The join of two graphs Hi and H2, denoted Hi + H2, is the graph formed by connecting 
each vertex of Hi to all of the vertices of H2 with edges. Then a (homogeneous) cyclic clan 
graph, denoted G[{Kr)m, jn] is a recursive graph composed of a set of m complete graphs 
Kr, such that the linkage between two adjacent pairs of Kr^s is a join (abbreviated jn). The 
Potts model partition function has the form ||5^ 



Z{G[{Kr)m,jn],q,v) = J2fid (W,r,dj)'" (6-1) 

d=0 j=l 



where /xq = 1 and 



where 



1{r) = Y[i(l - (6-3) 
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is the falling factorial in combinatorics. The numbers nz,cG{f,d) were determined in ||59[ . 

Let us apply the relation Z{G, g = 0, f ) = to this family. All of the Hd coefficients for 
d > 2 contain g as a factor and hence vanish identically for q = 0. Using /ii = g — 1, the 
condition Z{G, q = 0,v) =0 reduces to 

E i^Z,CG,r,0,r = E i^Z,CG,r,l,r . (6.4) 

i=i j=i 

Since this must hold for arbitrary m and since, in general, the relevant A's are nonvanishing, it 
follows that the cancellation must occur in a pairwise manner, and hence that nz,cG{f, 0) = 
nz,cG{f-,^)- This is in agreement with a result that we already derived in [^, viz., that 
^z,CG{fi 0) = nz^cGi^i 0) = 2''"^. So here one does not obtain any new sum rule. A similar 
comment applies for the chromatic polynomial P{G[{Kr)m) jn]., q), where np^cG{fj d) = 1 for 
all d (with < d < r) [0. 

7 Determination of Xp^Ly,d=Ly-i for Families of Lattice 
Strip Graphs 

7.1 Cyclic and Mobius Strips of the Square Lattice 

We have succeeded in determining the terms \p,Ly,d=Ly-i with coefficient c^'^y'^^ in the 
chromatic polynomial for the cyclic and Mobius strips of the square lattice. There are 
np{Ly,Ly — 1) = Ly of thcsc terms, by eq. (|2.2.22|) . Since all of the results in this section 
and the rest of the paper (including the appendix) refer to the chromatic polynomial, we 
shall use an abbreviated notation omitting the subscript P in the terms \p^Ly,d- For our 



calculation, we use the sieve method of to calculate the Ly x Ly transfer matrix for 



these A's with coefficient c*-'^""^-'. This transfer matrix is denoted T^^y''^\ To avoid awkward 
notation, we leave the type of lattice implicit. In subsequent sections, when the same symbol 
is used for strips of other lattices, it will be understood that its meaning is specific to those 
sections. It is convenient to extract a prefactor and define a reduced transfer matrix: 

Theorem 1 For the cyclic or Mobius strip of the square lattice with arbitrary width 
Ly and arbitrary length L^ = m, the eigenvalues XLy,d,j of for coefficient degree d = Ly — 1 
are 

K,'i=Ly-lJ = (-l)^'' + '(g - a,,,Ly,j) ,l<J<Ly (7.1.2) 
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where 



^sq,Ly,j = 1+4 COS^ 



{Ly + l-j)n 
2L„ 



, l<J<Ly. 



(7.1.3) 



Proof 

We begin with the following lemma: 

Lemma 1 For the cyclic strip of the square lattice with arbitrary width Ly and length 
= m, denoted sq[Ly x L^, , eye], the matrix T^^i^'^^'^^ for Lj^ > 2 is given by 



rp{Ly,Ly — l) 
^ 11 



(Ly,Ly — l) 



q-2 



(7.1.4) 



rp{Ly,Ly — l) 

33 



q-3 for 2<j<Ly-l 



(7.1.5) 



rp{Ly,Ly — l) 

^3 i+1 



rp{Ly,Ly — l) 

3 



-1 for l<j<Ly-l 
with other elements equal to zero. For Ly = 1, T*^^'°) = Ai^o = Q ~ ^■ 
Thus, for example, for Ly = 5, 


















f q-2 -1 

-1 q-3 -1 

-1 q-3 -1 
0-1 q-3 -1 
\ -I q-2 J 



(7.1.6) 



(7.1.7) 



We note that T^^^^y is a special case of a Jacobi matrix, in the terminology of linear 
algebra (where a Jacobi matrix A is defined as a matrix with the property that Aij = if 
\t~j\>2) 0. 

Proof One can construct the cyclic strip of the square lattice via the repetition of a 
subgraph which is a transverse slice consisting of the path graph with Ly vertices such that 
the linkage L between two successive such path graphs is the identity linkage, so that vertex 
1 of the first path graph is connected by an edge to vertex 1 of the next path graph, and 
so forth for all Ly vertices. We denote this linkage as L = {(1, 1), (2, 2), {Ly, Ly)}. Now 
applying the sieve method of we consider the functions a : V {1,2, that are 
proper g-coloring of the path (or tree) graph with Ly vertices, so that V is the vertex set 
V = {1,2, Ly}. Let us choose the basis so that 



N(/5) 



1 if a = (3 
otherwise 



(7.1.8) 
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Define the compatibility matrix as follows, 

/rp\ _ J 1 if is compatible with L (7 ^ q\ 

I- otherwise 

Denote [l\h] as the function that takes the value 1 for the coloring that assign color h to 
vertex 1, and takes the value otherwise. Similarly let [1, 2, ...\hi, h2, ...] refer to the coloring 
choice such that color hi is assigned to vertex 1, color /12 is assigned to vertex 2 and so forth. 
The coloring without restriction is denoted u. We have 

xcv 

where X refers to all subsets of the Ly vertices of the path graph, and |X| is the number of 
vertices in X. For example, for Ly = 3, we have 

T[a] = u-[l\ai] - [2\a2] - [3|a3] + [1, 2\ai, 02] + [1, 3|ai, a^] + [2, 3\a2, 0:3] - [1, 2, 3|ai, 0:2, as] . 

(7.1.11) 

Concentrating now on the case ci = — 1, we observe that the invariant subspace of the 
matrix T is spanned by 

[l,2,...,j - l,j + 1,.. .,1.3^1/11, /i2,...,/ij_i,/ij+i,...,/ilJ, l<j<Ly. (7.1.12) 

Therefore, 

T[l, 2, J - 1, J + 1, Ly\h, h2, hj.i, hL,] = T ( j = J2T[a] (7.1.13) 

\ H J H 

where H means ai = hi, 0^2 = /i2, Oj-i = ^j-i, ctj+i = ^j+i, ctLy = hiy, and 
hi, h2, hj-i, hj+i, hiy are different from each other. 

To obtain the matrix for level Ly — 1, we only have to consider the terms in the summation 
of eq. ( |7.1.iq ) with |X| = - 1 and |X| = L^. Now the first row of T^Ly^y-^) is obtained 
by applying T on [1, 2, Ly — l\hi, h2, /^Lj^-i]. The nonzero elements are 

{-l)^y-\{q - 1)[1, 2, Ly -l\hi, h2, hLy-i] + (-[1, 2, Ly -2,Ly\hi, /i2, /il,-2, /^l 

+ (-l)^-[l,2,...,L,-l|/li,/l2,...,/lL,-l] 

= (-l)^--^((g - 2)[1, 2, Ly -l\hi, h2, hL,-i] 

- [l,2,...,Ly-2,Ly\hi,h2,...,hLy-2,hLy]) (7.1. 
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and we have the corresponding result for the last row. Applying T on [1,2, — l,j + 
1, Ly\hi, /i2, hj^i, hj^i, h^y] for 2 < j < Ly — 1, we find that the nonzero elements 
are 

{-l)^y-\-[l, 2, J, J + 2, Ly\hi, /i2, hj, hL,] 
+ {q - 2)[1, 2, j - 1, j + 1, Ly\hi, /i2, /ij-i, /ij+i, hLy] 
+ (-[1, 2, j - 2, j, Lj^|/ii, /i2, hj_2, hj, /lij)) 

+ (-l)^^[l,2, - l,j + l,...,Ly|/li,/l2,---,/ij-l,/ii+l,---,/iLj 

= 2, J, J + 2, Ly\h, /i2, /ij, hL^] 

+ {q - 3)[1, 2, j - 1, J + 1, Ly\hi, /i2, /^i-i, /ij+i, hiy] 

-[1, 2, j - 2, j, /i2, /ij-2, /ij, /iij) (7.1.15) 

and, taking into account the relation ( [7.1.1|) , the lemma follows for the case of the cyclic 
strip. It was proved in that the A's for the cyclic and Mobius strips of the square lattice 
are identical to each other, and similarly for the triangular lattice. Finally, for the Ly = 1 
case, = g - 1 is obtained by an exphcit elementary calculation. This completes 

the proof of the lemma. □ 

Next, we have 

Lemma 2 The eigenvalues of the (reduced) transfer matrix T^^y'^y~^^ are q — cbLy,j for 

1 < J < Ly. 

Proof For the case Ly = 1, an elementary explicit calculation yields the result. For 
Ly > 2, consider the Ly x Ly matrix A{Ly) given by 

A{Ly)j ,+, = A{Ly)j+^, = -l for 1 < j < Ly - 1 (7.1.16) 

with all other elements equal to zero, and A{Ly) = for Ly = 1. Denote the characteristic 
polynomial of A{Ly) in terms of variable u as C[A(Ly)] = det(co'/ — A{Ly)), where / is the 
Ly X Ly identity matrix. This polynomial satisfies the recursion relation (for Ly > 2) 

C[A{Ly)] = uC[A{Ly - 1)] - C[A{Ly - 2)] (7.1.17) 
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where, for Ly = 2 we formally set C[A(0)] = 1. As in our earlier related work we 
observe that this is the same as the recursion relation satisfied by the Chebyshev polynomial 
of the second kind, Un{x), for x = uj/2, i.e., 

f/„(x) = 2xUn-lix) - Un-2{X) • (7.1.18) 

Solving the recursion relation ( 7.1.17 ), we obtain 

[— ] 

C[A{Ly)] = U,^{u/2) = X](-l)^ (^^ r^)(cu)^-^^- (7.1.19) 
where we again use the notation [u] to denote the integral part of u. Using the relation 



UJcos^) = !if*i±M (7.1.20) 

smd) 



we express this equivalently as 



= sinp + l)arccos(u;/2)] ^1) 
sin[arccos(u;/2)J 

It follows that the roots of the charactistic polynomial C[A{Ly)], i.e., the eigenvalues of the 
matrix A{Ly), are 

Uj = 2cos( ^ for l<j<Ly. (7.1.22) 

Now the reduced transfer matrix T^^y^y^^ can be written for Lj^ > 2 as 

f{Ly,Ly-i) = (^q_ 3)j _^ f(Ly,Ly-i) (7.1.23) 

and the matrix T^^y^i'"^^ is given by 

rp(Ly,Ly — l) rp{Ly,Ly—l) ^ (7 \ 2 

11 Ly Ly \ ' ' / 

fiLyJ^y-^) ^ f^^y^jy-^^ = _i for 1 < j < Ly - 1 (7.1.25) 

with all other elements 0. Let us denote the characteristic polynomial of T^^y'^y~^^ in terms 
of the variable u as C[T^^y''"y~^'']. We find that this characteristic polynomial satisfies the 
following recursion relation 

(j^f{Ly,Ly-i)^ = (cj - 2)C[A{Ly " l)] . (7.1.26) 
Using eq. ( |7.1.22| ), we find that the roots of C[T'^^y'^y~-^^] are given by 

coj = 2 cos(y^) for < j < Lj, - 1 . (7.1.27) 
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Combining eqs. (|7.1.23| ) and ( [7.1.27|) (and relabelling the roots so that the aLy,j increase 



monotonically as a function of j) then yields the result in the lemma. □ 



Finally, we use the result from |T^, |2^ that the Xp^a/^ for the cyclic and Mobius strips 
graphs G of a given lattice and width Ly are the same. The theorem then follows. □ 



We remark that eqs. (5.2) and (5.3) in are in accord with this structure of the 
transfer matrix. 

Several corollaries of Theorem 1 and the associated lemmas are of interest. Each of these 
applies to the cyclic or Mobius strips of the square lattice; we leave this implicit in the 
notation. Since the proofs are straightforward, we omit them. 

Corollary 1 The trace and determinant of T^^y'^y~^^ are 

T^j.(^fiLy,Ly-i)^ = 2 + Ly{q - 3) (7.1.28) 

and 

Ly 

det(f (^-^^-1)) = l[{q - a,,,Lyj) (7.1.29) 

where the asq,Ly,j were given in ( |7.1.3| ). 
Corollary 2 

For arbitrary Ly, g — 1 is an eigenvalue of T^^y^^^^^ with multiplicity 1, given by 

XLy,=i = q-l. (7.1.30) 

Corollary 3 

If = mod 2, i.e., Ly > 2 is even, then g — 3 is an eigenvalue of T^^f^""^^ with 
multiplicity 1 given by 

Al„j=i+(l,/2) = g - 3 . (7.1.31) 

Corollary 4 

If = mod 3, then q — 2 and g — 4 are eigenvalues of T^^y'^y~^\ each with multiplicity 
1, given by 

AL,.i=i+(L,/3) = g-2 (7.1.32) 
Al,j=i+(2L,/3) = g - 4 . (7.1.33) 

Corollary 5 

The roots asq,Ly,j monotonically increases from 1, for j = 1, to 1 + 4 cos^(7r/(2Lj^)) for 
j = Ly. As Ly — > oo, these roots become dense on the interval [1, 5]. 
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7.2 Self-Dual Cyclic Strips of the Square Lattice 

Theorem 2 For the self-dual cyclic strip of the square lattice with arbitrarily width Ly 
and length — m, the eigenvalues Xp^Ly,d,j for coefficient degree d — Ly — 1 are 

K4=Ly-i,j = - aG,,L,,i) ,l<j<Ly (7.2.1) 

where 

«G.,L„, = 1 + 4 cos^ .^<J<Ly. (7.2.2) 

Proof We begin with the following lemma: 

Lemma 3 For the self-dual cyclic strip of the square lattice with arbitrarily width Ly 
and length Lx = m, the reduced transfer matrix T^^^^^'^^ for > 2 is given by 

j,{Ly,Ly~i) ^ ^ _ 2 (7 2.3) 



T, 



{Ly,Ly — l) 



JJ 



q — 3 for 2 < j < Ly 



(7.2.4) 



r^{LyJ.y 1) ^ rp{Ly,Ly ^) foj. 1 < j < Ly - 1 (7.2.5) 

with other elements equal to zero. For Ly — 1, T^^'^^ = Ai,o = g — 2. 

Proof This proceeds in a manner similar to the proof given above, with the difference 
that the coloring is further constrained by the feature that all of the vertices on the upper 
side of the strip arc connected by edges to a single external vertex. This difference has the 
effect of replacing the g — 2 by g — 3 in the entry T^^^^ . 



□ 



An example for L^^ = 5 is 



f (5,4) 



/g-2 -1 

1 g- 3 
-1 



v 
















-1 

g-3 
-1 








-1 

g-3 
-1 









g-3 y 



Next, 

Lemma 4 



The eigenvalues of T^^^^y are given by 
q - aGn,Ly,j for l<j<Ly. 



(7.2.6) 



(7.2.7) 
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Proof For the case Ly = 1, an explicit calculation yields the result. For Ly > 2, 



consider the Ly x Ly matrix A{Ly) given by 



rp{Ly,Ly-l) _ _ _j_ J^iLy 



-1) 



where again I is the Ly x Ly identity matrix, and the matrix T^^y^^ is given by 



{Ly,Ly — l) 



11 



(7.2.8) 



(7.2.9) 



{Ly,Ly — l) 



j i+1 



{Ly,Ly — l) 



i+1 i 



-1 for l<3<Ly-l 



(7.2.10) 



with all other elements 0. Let us denote the characteristic polynomial of T^^y^y in terms 
of the variable to as C[T^^y'^y~^^]. We observe that this characteristic polynomial satisfies 



the recursion relations for Ly > 2 



{uj - l)C[A{Ly - 1)] - C[A{Ly - 2)] 



C[A{Ly)] - C[A{Ly - 1)] 



(7.2.11) 



where the matrix A{Ly) and its characteristic polynomial C[y4(Ly)] were given in Lemma 2, 
and we continue to use the formal definition C[A(0)] = 1. Using eq. ( |7.1.21| ), we find that 
(7[j'(-f^y.-C'y-i)] can be written as 

(Jlf{Ly,Ly-l)^ ^ sin[(L^ + l)arccos(a;/2)] _ sin[L^arccos(^/2)] 

sin[arccos(co'/2)] sin[arccos(co'/2)] 
2cos[(Lj^ + l/2)arccos(c<j/2)] sin[arccos(co'/2)/2] 

sin[arccos(co'/2)] 
cos[(Ly + l/2)arccos(c<j/2)] 
cos [arccos {u!/2)/2] 

Therefore, the roots of the characteristic polynomial ClT^^y^i^'^^] are given by 

■(2j-l)7r^ 



2 cos 



2L„ + 1 



for I < j < Ly . 



(7.2.12) 



(7.2.13) 



Combining eqs. ( [7.2.8| ) and ( |7.2.13| ) then yields the result in the lemma. □ 
Finally, combining Lemmas 3 and 4 yields the theorem. □ 

As is the case with asq^Lyj, the roots aGD,Ly,j become dense on the interval [1,5] as 
Ly ^ oo. It is also straightforward to derive the following corollary: 

Corollary 6 If = 1 mod 3, then g — 2 is a root of the reduced transfer matrix 
'f{Ly,Ly-i) self-dual cyclic strip of the square lattice of width Ly. 



28 



7.3 Cyclic and Mobius Strips of the Triangular Lattice 

We have also determined the corresponding Ly x Ly transfer matrix for the triangular lattice. 
As noted before, to avoid awkward notation, we shall use the same symbol for the transfer 
matrix T^^y'^\ but understand that it is different for this lattice strip than for the previous 
strips. We have the following theorem for T^^y'^>'~^\ 

Theorem 3 The general matrix T^^^'^^'^^ for the cyclic strip of the triangular lattice 
with > 2 is given by 

f(Ly,Ly-l)^^_^ (7.3.1) 

^{Ly,Ly-i) _ ^ _ 4 2 < z < - 1, l<j<i (7.3.2) 

fi!;'f'''^ =q-2 for l<j<Ly (7.3.3) 
f{LyM>-^) _ _^ 1 < J < Ly - 1 . (7.3.4) 

Proof We use the sieve formula again. Consider the strip of the triangular lattice of 
width Ly and arbitrary length to be constructed by starting with the corresponding cyclic 
strip of the square lattice and adding edges connecting the upper right and lower left vertices 
of the squares on the strip, so that the edge set is L = {(1, 1), (2, 1), (2, 2), (3, 2), ...(Ly, Ly)}. 
In this case, there is more than one edge connecting a vertex on one path graph forming 
a transverse slice to the adjacent path graph, so eq. ( [7.1. 101 ) has to be written in a more 
general form, 

TN= $^(-l)l^l$^[X|a.] (7.3.5) 

where we sum over all different sets of edges connecting the vertices of the right-hand path 
graph, X, to the vertices of the left-hand path graph for each such pair. As an example, for 
the Ly = 3 case, this is 

T[a] = u-[l\ai] - [l|a2] - Plas] - Pl^s] - [Sl^s] + [l,2|ai,a2] + [l,2|ai,a3] + [I,2|a2,a3] 

-F[1, 3|a;i, as] + [1, 3|a2, as] + [2, 3|a2, as] - [1, 2, 3|ai, a2, as] . (7.3.6) 

For (i = Ly — 1, the invariant subspace of T is again spanned by the terms given in eq. 
( [7.1.12| ). The nonzero elements in the first row of T^^y'^y~^^ are 



{-l)^y-\{q - 1)[1, 2, Ly - h2, /iL„_i] + (-[1, 2, Ly - h, hLy-i]) 
+ (-[1, 2, Ly - 2, Ly\h, h2, hLy-2, hLy])) + i-l f' [1, 2, Ly - l\hi, h2, hLy-i] 
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= {-lf^-\{q - 3)[1, 2, Ly - l\hi, /i2, hL,-i] - [1, 2, - /i2, /il,-2, /ilJ) ■ 

(7.3.7) 

For row j from 2 to L^^ — 1, the nonzero elements are 

2, - h2, K-i] + (g - 2)[1, 2, - h2, /iL,-i] 

+ (-[1, 2, Ly - l\hi, h2, hL,-i]) + ... 

+ (-[1, 2, J, J + 2, Ly\hi, h2, /ij, /ii+2, hiy]) 

+ {q - 2)[1, 2, j, j + 2, Ly\hi, /i2, /ij, /ii+2, /^lJ 

+ (-[1, 2, J, J + 2, Lj/|/ii, /t2, hj, hj+2, /iLj) 

+(q - 2)[1, 2, j + 1, /i2, hj-i, hj+i, /iij 

+(-[1, 2, J - 1, J + 1, Ly\hi, h2, /ij-i, hj+i, /ilJ) 

+(-[1, 2, - 2,j, Lj/|/ii, h2, /ij-2, /ij, /iij)) 

2, j - 1, J + 1, Ly\hi, /i2, -, Vi' ^ij 

= - 4)[1, 2, Ly - l\h^, h2, hL,-i] + ... 

+ (g - 4)[1, 2, j, j + 2, Lj^|/ii, h2, /ij, /ij+2, iiLy] 

+ (g-4)[l,2,...,j - l,j + l,...,Ly\hi,h2,...,hj_i,hj+i,...,hLy] 

-[1,2,...,J -2,j, ...,Lj^|/ii,/i2,...,/ij-2,^i,--,/iLj) ■ (7.3.8) 
The nonzero elements for the last row are 

2, Ly -l\hi, h2, /iL,-i] + {q- 1)[1, 2, -l\h,, h2, /iL,-i] + ... 

+ (-[1, 3, Ly\hi, hs, /ilJ) + (g - 1)[1, 3, Ly\hi, /ig, /i^J 
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+ (g - 1)[2, 3, Ly\h2, h, hL,]) + (-1)''^ [2, 3, Ly\h, 



l)'^--\{q - 2)[1, 2, ...,Ly-l\h, h2, + ... + (g - 2)[2, 3, L,|/i2, /is, /^lJ 



(7.3.9) 



The theorem then follows. □ 
Thus, for example, for Ly 



f (5,4) 



/g-3 -1 0\ 

g_4 g-4 -10 

g-4g-4g-4 -1 

g-4 g-4 g-4 g-4 -1 

\g-2 q-2 q-2 q-2 q-2 J 



(7.3.10) 



We also recall the result from [|1^, Q that the Ap^G.j's for the cyclic and Mobius strips graphs 
G of a given lattice and width Ly are the same. 

Here the eigenvalues of the cyclic triangular lattice do not have a simple linear form 
q '~'^Lyj as was the case for the eigenvalues for this degree c? = — 1 in the case of the cyclic 
and Mobius strips of the square lattice given in eqs. ( |7.1.2|) and ( |7.1.3|) or the eigenvalues 
for the self-dual strip of the square lattice given in eq. ( [7.2. 1|) with (|7.2.2D . 

However, we do find the following corollary, which is easily derived from the theorem. 

Corollary 7 For the cyclic strip of the triangular lattice 

det(f (^"'^"-i)) = (g - 2f{q - 3)^"-^ (7.3.11) 
^^{f{Ly,Ly-i)^^ ^ 3 ^ _ _ (7.3.12) 
For example, for Ly = 2, we have, for the eigenvalue of the reduced transfer matrix. 



2q-5± ^9-Aq 



for J = 1,2 



(7.3.13) 



in agreement with eqs. (5.12) and (5.13) of (taking account of eq. ( [7.1. 1|) ). For Ly = 3, 



we have 



A 



3,2,j 



^3,2,1 = q-2 

2q-7± ^25 - 8q ] for j = 2, 3 



(7.3.14) 
(7.3.15) 



in agreement our eqs. (2.16) and (2.17) in |[4^ . 
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8 Chromatic Polynomial for the Cyclic Strip of the 
Square Lattice of Width Ly = b 

8.1 General 

In this section we give our solution for the chromatic polynomial of the Ly x Lx cyclic strip 
of the square lattice with width Ly = 5. For beyond the first few degenerate cases, this 
has n = LxLy vertices and e = Lx{2Ly — 1) edges. The chromatic number is x = 2 for 
even and x = 3 for odd, independent of Ly. 

We list below the specific c^'^^^s that will appear in our results. 

c(^) = (g- I)(g3-6g2 + 9g_l) , (8.1.3) 

and 



For this family, using the general structural formulas derived in |^6[, we know that 

np(5,0) = M4 = 9 , np(5,l) =np(5,2) = 21 , np(5,3) = 13, np(5,4) = 5 (8.1.5) 

as well as np(5,5) = 1. Summing these or equivalently evaluating the general formula 
( p.2.24|) for the case Ly = 5, one has, for the total number of As^^^'s, Np^^^x = 70. We 
have calculated the chromatic polynomials by first computing the transfer matrices T^y^d for 
0<d<Ly. 
We find 

9 21 21 



P{sq[5xm,cyc.],q) = c^") ^(As.o,,)"^ + c« ^^(As,!,,)"^ + c^^) 5^(A5,2 ^™ 



j=i j=i j=i 

13 5 

+ c^''> E(^5,3,,)'" + c(^) 5^(A5,4,,)'" + c^'\X5,,r (8.1.6) 

where A5 o,j for j = 1, 2 and 3 < j < 9 are roots of equations of degree 2 and 7 in q, X^^ij 
for 1 < j < 8 and 9 < j < 21 are roots of equations of respective degrees 8 and 13, A5,2,j 
for 1 < j < 9 and 10 < j < 21 are roots of equations of respective degrees 9 and 12, and 
As.sj- for 1 < j < 5 and 6 < j < 13 are roots of equations of respective degrees 5 and 8. We 
discuss these A's next. 
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8.2 A5,oj 

The quadratic equations for factorize over the field Q[\/5], so that these eigenvalues 
are polynomials in q and the elements of this field: 



A 



5,0,j 



2q^ - 13q^ + 28q - 19 ± (g - l)(g - 3)^5 



for J = 1,2 



^.2.1) 



The Asoj, 3 < j < 9 are identical to the A's for the Ly = 5 strip with free boundary 
conditions, and the degree-7 equation for these is 

7 

+ Y.^s,,ffA'~' = (8.2.2) 



where the expressions bsq5FF,j were given in eqs. (A.20)-(A.26) of [p]) . 



8.3 ^5,d,j, — 1, 2 

The with d = 1 and d = 2 and 1 < j < 21 are too complicated to present here. The 
computer files are available on request from the authors. 



8.4 A5,3,j- 

The A. 3j- are the eigenvalues of the transfer matrix 



where 
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.4.1) 




til = 


-g^ + 4g - 5, 


^22 


= -(g 


-2)(g 


-3), 


^33 = 


-g^ + 5g - 8 


(8 


.4.2) 





t44 = -(g-2)(g-3), t55 = -(g-3)2, tm = -iq-2f 
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i77 = -(g-2)(g-3), t88 = -g' + 5g - 8, tgg = -(g - 2)(g - 3) 
^1010 = -g^ + 4g - 5 . 



As noted above, the characteristic polynomial of this matrix factorizes into degree-8 and 
degree-5 polynomials. 



8.5 A5,4,j and A5,5 

As a special case of our general formula ( [7.1.2| ) with ( [7.1.3|) for XLy,Ly-i,j, of which there are 
Ly eigenvalues, we have 

^5,4,1 = g ~ 0.sq,5,l = g — 1 



5 — y/b 

A5,4,2 = g - asg,5,2 = g " | ] = Q - 1-38196. 



^5,4,3 = g - asg,5,3 = Q - ^ ^ j = Q- 2.38196... 

A5,4,4 = g - asq,5A = Q - ^~~y~~^ = ^ " 3.61803... 
A5,4,5 = Q- «5g,5,5 = 9 - I ^"^y^ 1 = ^ ~ 4.61803... 



8.5.1) 
8.5.2) 

8.5.3) 

8.5.4) 

8.5.5) 



We note the connection between the appearance of the algebraic numbers a + by/5, i.e., 
elements of the field Q[v^] in both the X^^aj for 2 < j < 5 and in X^^/ for £=1,2. Indeed, 
the term involving \/5 in As^o.i and A5^o,2 is proportional to (g — l)(g — 3), which is A4^o,i- 
Our general formula ( p.2.25| ) yields 

A5,5 = -1 • (8.5.6) 



9 Chromatic Polynomial for the Mobius Strip of the 
Square Lattice of Width Ly = b 

We obtain this by applying the transformation rules derived in to go from the chromatic 
polynomial for the strip with width Ly and arbitrary length with {FBCy, PBC^) = cyclic 
boundary conditions to the corresponding strip with {FBCy, TPBC^) = Mobius boundary 
conditions. For the chromatic polynomial of the Mobius strip with width L^, = 5, we find 

7 9 2 12 

P(.g[5 X m,M6],g) = [ 5^(A5,o,+2)"^ + ^(^5,2,,)™ - E^^^a.)'" " E(^5,2,+9)™ 

j=i j=i j=i j=i 
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Expressing this in the general form ( 2.2.31| ), we have 



16 14 15 

P{sq[5xm,Mb],q) = c^") [X](A5,o,+,,)™ - ^^(Asa-,.)"^] + c^'^ [ E^^^ah 

i=i i=i 

i=i j=i 



where 



+ c(^HA5,5) 



A5,0H 



A. 



J - -^5,o,i+2 for 1 < j < 7 



A5,0, 



A, 



-7 for 8 < j < 16 



A5,o,-,i = A5,o,j for 1 < J < 2 



A5,0,- 



A 



^7 for 3 < J < 14 



A5,i,+,i = A5,ij+8 for 1 < j < 13 

A5,l,+,14 = ^5A,2 , ^5,1,+,15 = ^5,4,4 

= fo^ 1 < j < 8 

A5,l,-,9 = A5,4_i, A5^i__^io = A5^4^3, As^i^^^n = A5_4^5 
A5,2,+j = A5,3j+5 for 1 < j < 8 

A5,2 -,i = A5,3j for 1 < J < 5 



A 



(9.1) 



(9.2) 



5,3,+J 



\5,5 



(9.3 
(9.4 
(9.5 
(9.6 
(9.7 
(9.8 
(9.9 
(9.10 
(9.11 
(9.12 
(9.13 



10 Locus B for Ly = 5 Cyclic/Mobius Strips of the 
Square Lattice 

In this section we present the singular boundary B across which W{q) is singular in the 
complex q plane for the oo limit of the cyclic strip of the square lattice with width 



Ly = 5. We also recall that that B is the same for cyclic and Mobius strips [jT6|, ^ ^ 



as was proved in pO|. A plot of B is shown in Fig. 
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Im(q)0 




Re(q) 



Figure 1: Locus B for the oo limit of the strip of the square lattice of width Ly = 5 strip with 

{FBCy,{T)PBCx) = cyclic (equiv. Mobius) boundary conditions. For comparison, chromatic zeros calcu- 
lated for the strip length = 20 (i.e., n = 100 vertices) are shown. 
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The locus B is comprised of closed curves that separate the q plane into several regions. 
This locus crosses the real axis at g = 0, g = 2, and a maximal value, Qc, which is 



ge = 2.582385... (10.1) 

The region Ri includes the real axis for q > qc and q < and extends outward to complex 
infinity from the outer envelope of B. The region i?2 includes the real segment 2 < q < q^, 
while region R3 includes the real segment < g < 2. In region Ri, the dominant A is the 



root of the degree-7 equation obtained from the appendix of with the largest magnitude, 
which we label Xm. In region R2, the dominant is the root of the degree-12 equation 
with the largest magnitude, which we label Xfi2- Hence, qc is the solution of the equation 
of the degeneracy \Xri\ = \Xr2\- The value of gc that we have obtained for the present 
strip may be compared with the values that we obtained previously for the cyclic (Mobius) 
Ly X 00 strips of the square lattice with smaller widths, namely, qc = 2 for Ly = 1,2 ||16|| , 



qc = 2.33654... for Ly = 3 gg, and qc = 2.492845... for Ly = A We calculate that 
W{sq[5 X 00, FBCy, {T)PBC^], q) = 1.317594 at g^, as given in eq. ( [lOlD . In region R3, the 
dominant Xj is the root of the degree-13 equation with the largest magnitude, which we label 
Aij3. There are also complex-conjugate regions centered at approximately g = 2.7 ± 0.8i; 
we denote these as Ri,R\. The dominant Xj in these regions is the root of the degree-13 
equation of maximal magnitude here; this is denoted Xji4^. We have 

W = (Am)'/' , for g G i?i (10.2) 

\W\ = \Xr,\^/\ for gei?, , j = 2,3,4. (10.3) 

(In regions other than only the magnitude \W\ can be determined unambiguously ||16[|.) 
Our previous results such as [^, showed that there can also be very small sliver regions 
in the complex g plane; we have not made an exhaustive search for these. As we found 
previously for = 3 [^, |36| and Ly = A [^, the locus B has support for Re{q) < as well 



as Re{q) > 0. 

Fig. |1| also shows a comparison of the chromatic zeros calculated for a long finite cyclic 
strip with L^ = 20, i.e., n = 100 vertices, versus the asymptotic locus B. One sees that, 
aside from the isolated real zero at g = 1, these chromatic zeros lie reasonably close to the 
locus B. On the curve forming B passing through g = and q = qc, the density is highest 
on the right-hand side and somewhat lower in the region of g = 0. The density of chromatic 
zeros on the part of B passing through g = 2 and ending at upper and lower triple points 
is very low, while one observes an intermediate density on the complex-conjugate pair of 
bulb-like curves protruding to the right. These features are qualitatively the same as we 
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found for Ly = 3 (Fig. 1 of |3^) and for Ly = 4 (Fig. 1 of ^Tl). The differences in densities 
on different portions of B were less pronounced for Ly = 2 (Fig. 1 of [0). The elementary 
case = 1 is special in two respects: (i) the complex chromatic zeros lie exactly B for finite 



as well as infinite L^., and (ii) their density is constant on i3 [16 

11 Discussion 

Here we make several further comments on these results for B 

1. As in our earlier work, we characterize a strip graph as containing a global circuit if 
it contains a cycle along the longitudinal direction whose length goes to infinity as the 
length L;j. — > oo; in practice, this is equivalent to the property that the lattice strip 
graph has periodic longitudinal boundary conditions. For all of the strips of the square 
lattice containing global circuits that we have studied, the locus B encloses regions of 
the q plane including certain intervals on the real axis and passes through g = and 
g = 2 as well as other possible points, depending on the family. Note that the presence 
of global circuits is a sufficient, but not necessary, condition for B to enclose regions, as 
was shown in |3l[] (see Fig. 4 of that work). Our present results for the square lattice 



are in accord with, and strengthen the evidence for, the inference (conjecture) [T^, [19 
that 

Bd {q = 0, 2} for sq[Ly,FBCy,{T)PBC^] V > 1 . (11.1) 

(Of course Ly = 1 graphs with {FBCy,TPBCx) and {FBCy, PBCx) boundary condi- 
tions are identical.) 

2. The crossing of B at the point g = 2 for the (infinite- length limit of) strips with global 
circuits nicely signals the existence of a zero-temperature critical point in the Ising 
antiferromagnet (equivalent to the Ising ferromagnet on bipartite graphs). This has 
been discussed in |T9| in the context of exact solutions for finite-temperature Potts 
model partition functions on the Ly = 2 cyclic and Mobius strips of the square lattice. 
In contrast, this connection is not, in general, present for strips with free longitudinal 
boundary conditions since B does not, in general, pass through q = 2. Furthermore, for 
the strips without global circuits, there is no indication of any motion of the respective 
loci B toward q = 2 a.s Ly increases. 

3. For cyclic strips, we note a correlation between the coefficient cq^j of the respective 
dominant Ag^j's in regions that include intervals of the real axis. Before, it was 
shown [|l^ that the cg^j of the dominant XgsJ region Ri including the real intervals 
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q > qc{{G}) and g < is c^^-' = 1, where the c*^*^^ were given in eqs. ( p.2.1|) . We observe 



further that the cg^j that multiphes the dominant XgsJ the region containing the 
intervals < g < 2 is c^^\ For the cychc Ly = 3 and Ly = A strips, there is also another 
region containing an interval 2 < q < on the real axis, where qc — 2.34 and 2.49 for 
Ly = 3,4; in this region, we find that the cg^j multiplying the dominant XgsJ is c*^^^. 
Our present results again agree with this conjecture. 

4. For the L^, — oo limit of all of the strips of the square lattice containing global 
circuits, a qc is defined, and our results for the cyclic and Mobius strips with widths 
from Ly = 1 through Ly = 4 indicate that qc is a non-decreasing function of Ly in 
these cases. The same behavior was found for the strips of the triangular lattice with 



Ly = 2 (and subsequently also Ly = 3,4 This motivated the inference 



(conjecture) that qc is a non-decreasing function of Ly for strips of regular lattices with 



{FBCy, {T)PBCx) boundary conditions |jT8[, and our present results strengthen the 
support for this inference. Given that, as Ly oo, qc reaches a limit, which is the qc 
for the 2D lattice of the specified type (square, triangular, etc.), this inference leads to 
the following inequality: 

g,(A, Ly X oo, BCy, {T)PBC.,) < qdA) . (11.2) 

This is a non-decreasing function of Ly. Our current results are in agreement with this 
conjecture. 

5. One interesting feature of our new results is that the outermost envelope of the curves 
on B for the Ly = 5 strip of the square lattice do not lie outside the envelope for B for 
the Ly = 4 strip. This is evident from the comparison shown in Fig. ^. For reference, 
we also show the analogous pairwise comparisons for widths Ly = 2,3 and Ly = 3,4. 
For each of these, it was true that as Ly increases, the outer envelope of the locus B 
moves outward in a manner such that B for width Ly + 1 encloses the locus B for Ly, 
allowing for equality at some points, in particular, g = |TB[ . This is also true for the 



pairwise comparison of the loci B for Ly = 1 and = 2, as is clear from Fig. ^ and 
the fact that the locus for = 1 is the circle |g — 1| = 1. However, our new results 
show that this behavior does not persist in the comparison oi Ly = 4 and Ly = 5 and 
hence is not general. 
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Figure 2: Comparison of loci B for the ^ oo limits of the strips of the square lattice of widths Ly = 4 
(light curves) and Ly = 5 (heavy curves) with {FBCy,{T)PBCx) = cyclic (equiv. Mobius) boundary 
conditions. 



40 



T 




Re(q) 



Figure 3: Comparison of loci B for the ^ oo limits of the strips of the square lattice of widths Ly = 2 
(light curves) and -Lj, = 3 (heavy curves) with {FBCy,{T)PBCx) = cyclic (equiv. Mobius) boundary 
conditions. 
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Figure 4: Comparison of loci B for the ^ oo limits of the strips of the square lattice of widths Ly = 3 
(light curves) and Ly = A (heavy curves) with {FBCy,{T)PBCx) = cyclic (equiv. Mobius) boundary 
conditions. 
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12 Appendix: Structure of Chromatic Chromatic Poly- 
nomials for Cyclic and Mobius Strips of the Square 
Lattice 

12.1 Ly = 2 Cyclic Strip 

We list some known results to illustrate our general notation, for comparison with our new 
calculation of the Ly = 5 case. The chromatic polynomial of the cyclic strip with width 
Ly = 2 is U 



P{2xm,cyc.,q) = cW(A2,o,i)'" + c« J](A2,i,,)"^ + c(2nA2,2) 



i=i 



iq' -3q + Sr + (g - 1) [(1 - qr + (3 - g)"^] + (q' - 3g + 1) 



where 



A2,o,i = - 3g + 3 
^2,1,1 = -{q - 0-2,1) = 1 - g 
A2,i,2 = -(g - 02,2) = 3 - g 

A2,2 = 1 • 



;i2.i.i) 

;i2.i.2) 
;i2.i.3) 

12.1.4) 

;i2.i.5) 



12.2 Ly = 2 Mobius Strip 



The chromatic polynomial of the Mobius (Mb) strip with width = 2 is p2 



P{2xm,Mb,q) = c^") [(Aa.ci)'" - (A2,2r] + c« ^^(-ly (A2,i,,r 

_ 3g + 3)- _ 1] + (g _ 1) J _ (1 _ g)- + (3 _ g)- (12.2.1) 



Expressing this in the general form (|2.2.31|), we have 



P(2 X m,Mb,q) = c^'^ {X2,o,+,T - (A2,o,-,i)^ 
where 



+ d 



(1) 



-^2,0,+,! — -^2,0,: 



(A2,i,+,i)"^ - (A2,i,-,i)'"J (12.2.2) 

(12.2.3) 
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-^2,0,-,! = ^2,2 
-^2,1, +,1 = -^2,1,2 
•^2,1,-,! = -^2,1,1 • 



12.3 Ly = 3 Cyclic Strip 



The chromatic polynomial for the cyclic strip of width L^^ = 3 is |3g, |36 



2 4 3 

P(3 X m, eye, q) = c(°) $^(A3,o,.)'" + c^'^ ^(^3,1,,)'" + c^'^ E(^3,2,,0 

j=i j=i j=i 



where X^^j for j = 1, 2 are the solutions to a quadratic equation given in PB|] , 

A3,i,i = -(g - 2)2 



A3,2,2 = q- a3,2 = g - 2 

^3,2,3 = q - (33,3 = g - 4 . 



Finally, 



A 



3,3 



-1 . 



12.4 Ly = 3 Mobius Strip 



The chromatic polynomial for the Mobius strip of width = 3 is 



Pi3xm,Mb,q) = c^") [ 5^(A3,o,r + $^(-1^ (A3,2, 

i=i j=i 

3 



+ C(l) [^(A3,l,, + l)™ - (A3,l,l)"^J + C^'\X3,3y 



(12.2.4) 
(12.2.5) 
(12.2.6) 



(12.3.1) 



(12.3.2) 



the A3 ij- for j = 4, 5, 6 are roots of a cubic equation given in and A3_2,j = Q — o.3,j, with 
the specific values 

A3,2,i = q- 03,1 = g - 1 



(12.3.3) 
(12.3.4) 
(12.3.5) 

(12.3.6) 



(12.4.1) 



Expressing this in the general form ( |2.2.31| ) 



3 2 
P(3 X m, Mb, q) = c^") [ 5^(A3,o,+,,)'" - $^(A3,o,+,,; 

i=i i=i 
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+ c^'^ - (A3,i,-,i)"^] + c^'\X3,3r (12.4.2) 



where 



A3,o,+,i = h,o,j for j 

^3,0,+,3 = ^3,2,2 
-^3,0,-, 1 = ^3,2,1 



1,2 



^3,0,-,2 — ^3,2,3 



A3,i,+,i = for j 

^3,1,-,! = A3, 1,1 
^3,2,+,! = ^3,3 . 



1,2,3 



12.5 Ly = 4 Cyclic Strip 



The chromatic polynomial of the cyclic strip of width L^^ = 4 is pi 

4 9 



P(4 X m, eye] 



4 



+ c^'^Y.^\,3,T+c^'\\,,r 



where 



(12.4.3) 
(12.4.4) 
(12.4.5) 
(12.4.6) 
(12.4.7) 
(12.4.8) 
(12.4.9) 



(12.5.1) 



A4,o,i = (?-l)(9-3) (12.5.2) 

A4,oj for i = 2, 3, 4 are roots of a cubic equation, XA,i,j for 1 < j < 4 and 5 < j < 9 are roots 
of equations of respective degrees 4 and 5, \i,2,j for 1 < j < 3 and 4 < j < 8 are roots of 



equations of respective degrees 3 and 5, all given in and A4,3j = — g, with specific 
values 

A4,3,i = -(<?- «4,i) = 1 - g (12.5.3) 

A4,3,2 = -(9-04,2) = 3- V2-g (12.5.4) 

A3,i,3 = -{q - (34,3) = 3 - g (12.5.5) 

A4,3,4 = -(5-04,4) = 3 + v^-g . (12.5.6) 

Finally, 

A4,4 = 1 . (12.5.7) 
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12.6 Ly = 4: Mobius Strip 

The chromatic polynomial for the Mobius strip of width L^^ = 4 is [^] 

3 3 5 

P{4xm,Mb) = c(°)[5^(A4,o,,+i)'" + $^(A4,2,,r-(A4,o,i)'"-5Z(A4,2„ 



2j+3j 



5 



4 



i=i i=i 



+ c(2)^(-iy(A4,3,,) 



(12.6.1) 



where 



the general form (2.2.31|), we have 




6 

P(4xm,M6) = c(°)[5^(A4,o,+,i)" 

i=i 


6 

- y^(A4,o-,j)" 


5 

j=i 


- E(A4,i,-,,)"^" 


+ C(^^[E(A4,2,+.)" 

i=i 


- E(A4,2,-,,)"^" 



A4,o,+j = A4,oj+i for 1 < j < 3 
A4,o,+j = A4,2,i--3 for 4 < j < 6 

A4,0,-,l = A4_o,l 



A4,0.-,7 ~ A, 



for 2 < j < 6 



A4,i,+j = A4,i,i+4 for 1 < j < 5 
A4,i -J = A4,i,j for 1 < j < 4 

A4,l,-,5 = A4^4 
A4,2,+,l = A4,3^2 , A4,2,+,2 = A4,3,4 
A4,2,-,l = A4,3^i , A4,2,-,2 = A4,3^3 . 



(12.6.2) 

(12.6.3) 
(12.6.4) 
(12.6.5) 
(12.6.6) 
(12.6.7) 
(12.6.8) 
(12.6.9) 
(12.6.10) 
(12.6.11) 
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